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Abstract 

The effect of unitary noise on the discrete one-dimensional quantum walk is studied using com- 
puter simulations. For the noiseless quantum walk, starting at the origin (n = 0) at time t = 0, the 
position distribution Pt(n) at time t is very different from the Gaussian distribution obtained for 
the classical random walk. Furthermore, its standard deviation, a it) scales as a{t) ~ t, unlike the 
classical random walk for which a(t) ~ \ft. It is shown that when the quantum walk is exposed 
to unitary noise, it exhibits a crossover from quantum behavior for short times to classical-like 
behavior for long times. The crossover time is found to be T ~ a~ 2 where a is the standard 
deviation of the noise. 

PACS numbers: PACS: 03.67.Lx,05.40.Fb 
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I. INTRODUCTION 

Random walk models describe a great variety of diffusion phenomena in physical systems. 
Such phenomena include the diffusion of particles in a fluid (Brownian motion), the motion 
of vacancies in a crystal and of atoms on a crystalline surface. Related models are also used 
to describe the spatial structure of systems such as polymer chains. 

A random walker on a lattice hops at each time step from its present site to one of 
its nearest neighbors. The hopping direction is picked randomly, with, for example, the 
probability to hop to each of c nearest neighbors given by 1/c. Random walk models on 
discrete lattices as well as in the continuum have been studied extensively. 

Consider a random walker on a one dimensional lattice. Denote the probability to find 
it at site n at time t by Pt(n). The time evolution of Pt(n) is described by the recursion 
equation 



For a walker starting from the origin at t = 0, the probabilities, Pt{n), of the walker to be 
at site n at time t, are given by the components of the Pascal triangle, namely 



and Pt{n) = for \n\ > t. Note that at even times only the even sites can be occupied, while 
at odd times only the odd sites can be occupied. A continuum description of the random 
walk can be obtained when the lattice constant An — > and the time step At — > 0, such 
that D = An 2 j (2 At) converges to a finite value. The parameter D is called the diffusion 
coefficient and it quantifies the rate in which the random walker moves. In this limit the 
random walk can be described by the diffusion equation 



where n and t are continuous variables. The solution of this equation, approached by Eq. 
(J2J) in this limit, is the Gaussian probability distribution 
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with D — 1/2. The standard deviation of this distribution, a(t) = \^2Dt, thus takes the 
form <r(t) = \/t. 

Recently, quantum analogues of the classical random walk model have been studied on 
the one dimensional lattice Q, 0, [i , ^ , 0] as well as on more general graphs . Consider 
the quantum walk on the one dimensional lattice. Unlike the classical random walker that 
occupies a single site at a time, the quantum walker can be in an extended state. This 
state is a superposition of all the basis states |n), n = 0, ±1, ±2, . . . , where |n) is the state 
in which the walker is located at site n. In addition to the spatial degree of freedom the 
walker has a chirality qubit, which can be in a superposition of the states \R) and \L), and 
determines the direction of the next hop. Each move of the walker consists of two unitary 
operations. The first one is a unitary transform, taken here to be the Hadamard transform, 
on the chirality qubit (expressed in the standard basis). The second one is the actual move, 
in which a walker with chirality \R) moves to the right and a walker with chirality \L) moves 
to the left. This definition resembles the classical walk in the sense that a walker in a basis 
state \n)\s) (s = R or L) at time t, will have equal probabilities to be found in the n + 1 
and n — 1 sites if a measurement is performed at time t + 1. However, the coherent motion 
of the state gives rise to constructive and destructive interference that strongly modifies the 
emerging probability distribution. As a result, a quantum walk starting from the origin at 
time t — moves away faster than the classical random walk. The standard deviation a(t) 
of the probability distribution Pt{n) for the quantum walk increases with time according to 
a{t) ~ t, while for the classical random walk a{t) ~ \fi. 

One of the main difficulties in the experimental implementation of quantum algorithms 
is the sensitivity of the quantum systems to noise [g, 0, and decoherence 
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problem of decoherence can be tackled by using quantum error correction methods 
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161 ] as well as decoherence free sub-spaces |l7|, [la iTjJ |2C| . However, these methods 
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require significant overhead, making the quantum circuits more complicated. It is thus 
useful to examine the effect of decoherence on various quantum algorithms, implemented 
on unprotected quantum circuits [2^. The quantum walk model is particularly suitable for 



this task since its behavior is qualitatively different from the classical random wa 
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effect of decoherence on the quantum walk has been studied recently 
It was found that the quantum walk is highly sensitive to decoherence, namely, even weak 
decoherence gives rise to classical-like behavior in the long time limit. Decoherence is a 
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fundamental problem because it is a result of the unavoidable interaction of the system with 
the environment, that makes the two of them entangled. This interaction can be described 
by a variety of noise models, bringing the system into a mixed state, by essentially tracing 
out the environment degrees of freedom. 

A related problem is the effect of unitary noise, which appears as a result of fluctuations 
and drifts in the generating Hamiltonian of any given unitary operation 27[. This noise 
tends to reduce the performance of quantum devices, but does not cause an entanglement 
with the environment, namely the system remains in a pure state. 

In this paper we analyze the effect of unitary noise on the discrete one-dimensional 
quantum walk model, using computer simulations. We find that even a tiny noise level 
will eventually induce a crossover from the quantum walk into a classical-like behavior, 
characterized by a(t) ~ \ft. The crossover time T is calculated numerically. It is found that 
T ~ a~ 2 , where a is the standard deviation of the noise distribution. 

The paper is organized as follows. The quantum walk model on the one-dimensional 
lattice is described in Sec. II. The noisy quantum walk is introduced in Sec. III. Simulations 
and results are presented in Sec. IV, followed by a summary in Sec. V. 

II. THE QUANTUM WALK MODEL 

Consider a quantum walker on a one-dimensional lattice. The state of the system at time 

t is 

oo R 
n=— oo s=L 

where the amplitudes a n>s are complex numbers that satisfy X^nJ a n,s| 2 = 1- The 
vector \n), where n = 0,±1,±2,..., represents the position of the walker on the lattice. 
The state vector \s) = \R) or \L) represents the chirality degree of freedom. The chirality 
consists of a single qubit, and its state determines the coefficients a niS (t + 1). Operations in 
the chirality space in the quantum model, replace the randomized decision that appears in 
the classical walk model and determines the hopping direction. When measuring the walker 
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position at a certain time t, the probability to find the walker at site n is given by: 



Pt(n) = k,;?(t)| 2 + \a n , L (t)\ 2 . (6) 
The time evolution of the quantum walk is expressed by 

m+i))=Q \mh m 

where \ip(t)) is the quantum state at time t. The quantum walker's step is defined by 

Qo = fU , (8) 

where U is a unitary operator that applies only in the chirality space and takes the role of 
the "coin" in the classical random walk. Here we focus the case of the Hadamard walk, in 
which U = I <S> wo, where 

is the Hadamard operator. The states of the chirality qubit are expressed in the standard 
basis 

|a>= ; \l)= ^ | , (io) 

and / is the spatial identity operator. The computational basis state vectors are transformed 
by the Hadamard operator according to 

U Q \n)\R) = -L\ n )(\R) + \L)) 

U \n)\L) = ±=\n){\R)-\L)). (11) 

The translation operator T then performs the walker's move according to the chirality state 
such that: 



f\n)\R) = \n + l)\R) 
f\n)\L) = |n-l)|L). 



(12) 



These operators yield the following recursion equations for the amplitudes of the quantum 

state 



a n , L (t + l) = -^=[a n+1>R (t) -a n+1>L (t)}. (13) 

For a given initial state, the recursion equations provide the probability of finding the walker 
at any site n at time t. Here we focus on the case that at t — the walker is located at the 
origin, namely, a n s (0) = for all n ^ (and s = L,R). 



III. THE NOISY QUANTUM WALK 

The quantum walk, like other quantum computing systems, may be affected by unitary 
noise. This noise is due to fluctuations and drifts in the parameters of the quantum Hamil- 
tonian of the system. The perturbed Hamiltonian is still Hermitian and therefore generates 
time evolution operators that are unitary, but now include a stochastic part. Formally, one 
can write a noisy unitary operator as 

U = U e iA , (14) 

where A is a stochastic Hermitian operator determined by the perturbation and U is the 
time evolution operator of the original quantum process without the perturbation. 

The Hadamard walk model consists of Hadamard gates that act in the chirality space at 
each time step of the walker. Here we consider the effect of unitary noise in the Hadamard 
operator on the quantum walker. To this end we describe a move of the noisy walker at 
time t as: 

m+i))=Q(t)\m)- (is) 

where 

Q(t) =f{I®w e l&{t) ). (16) 
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The operator a(t) is a stochastic and Hermitian operator that applies in the single-qubit 
chirality space at time t. It can be expanded in the basis of the Pauli operators such that 



a(t) = ai(t)<7i + a 2 (t)a 2 + a 3 (t)a 3 , (17) 

where a k (t), k = 1,2,3 are real stochastic variables, and the cr's are the Pauli operators, 
represented in the spin basis as: 

*='!!)' * 2= {°l)' *° = {l-i l - (18) 

Although, in principle, the 2x2 identity operator should also be included in the expansion 
(|17|) . it is omitted because it can only change the overall phase in the quantum walk process, 
and does not have any effect on the walker's measurement probabilities. 

We will now analyze the properties of the noisy Hadamard walk according to the charac- 
teristics of the real stochastic variables a k (t), k = 1, 2 and 3. In the analysis we assume that 
there is no correlation between different noise components k and k' ^ k as well as between 
different walker steps at times t and t' ^t, namely 

(a k (t)a k ,(t')) = 8 k , k >8 t ,t'Ot 2 , (19) 

where 5 kjk i is the Kronecker delta function. Furthermore, we focus on the isotropic case in 
which «i, «2 and are taken from the same distribution p(a), which is unbiased, namely 

(a fc (t))=0, = 1,2,3, (20) 

with standard deviation a. 

To analyze the time evolution of the Hadamard walk in the presence of unbiased and 
isotropic unitary noise, we have performed computer simulations of the model. In these 
simulations the Hadamard gate is replaced by a noisy one, so the walker step at time t is 
given by Eq. (|16)). The Pauli coefficients a\(t), atzit) and a$(t) are taken from a Gaussian 
distribution with zero average and a standard deviation a. The magnitude of the isotropic 
noise is characterized by the values of a. To obtain proper statistics the results were averaged 
over a sufficient number of runs. For each value of a, we have applied 10, 000 steps of the 
noisy Hadamard walk for at least 200 runs in the case of weak noise (a < 0.07) and more 
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than 4,000 runs in case of strong noise (a > 0.07). The distribution Pt{n) of the walker's 
position was then examined and its moments are calculated. 

The probability to measure the noisy Hadamard walker in site n at time t is given by: 

(Pt(n)) a = (\a n , R (t)\ 2 ) a + (|a n , L (t)| 2 ) a , (21) 

where a nj n(t) and a„,L(t) are amplitudes of the walker's state \ip(t)) [see Eq. ©]. Here () Q 
denotes the averaging over the noise, taken from a Gaussian distribution with a zero average 
and standard deviation a. 

Here we consider only initial conditions for which at t = 0, the walker is located at the 
origin, namely a njS (0) = for all n ^ 0. For these initial states it is guaranteed that Pt{n) = 
for all \n\ > t. Therefore the first and the second moments of the spatial distribution at 
time t (which are both averaged over the noise) are given by: 

t 

Wfi) a = n{Pt{n)) a (22) 

n=—t 

and 

t 

W))a = n2 ( p t(n))a- (23) 

n=— t 

The standard deviation a a (t) of the spatial distribution of the noisy Hadamard walk at time 
t is given by 

= - (W)) a 2 - (24) 



IV. SIMULATIONS AND RESULTS 



To analyze the effect of noise on the Hadamard walk we have performed direct computer 
simulations of the system, using Eq. (|16p. The initial condition used in the simulations is 

= ^(|0>|i?) + z|0)|L)). (25) 

For this state the probability distribution Pt(n) of the Hadamard walk without noise turns 
out to be symmetric around the origin |5|, namely Pt(—n) = Pt{n). In Fig. ^we present the 
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probability distribution P t (n), n = -t, -t + 2, . . . , t, for t = 250 (a), 1000 (b), and 10, 000 
(c) steps. The distribution exhibits two main peaks, on opposite sides of the origin, that 
move away from each other as time evolves. In the central part it forms a plateau, in which 
Pt{ n ) ~ It also exhibits wild oscillations under the envelope, that decays towards the 
origin and rises towards the two peaks. This distribution is clearly very different from the 
Gaussian distribution obtained for the classical random walk. The formation of Pt{n) can 
be understood as a result of constructive interference in the front of the propagating wave 
function, while the center is dominated by destructive interference. 

When unitary noise is added, it perturbs the structure produced by constructive and 
destructive interference. The distribution Pt{n) for a weak noise level (a = 0.025) is shown 
in Fig. H for t = 250 (a), 1000 (b), and 10,000 (c) steps. The effects of the competition 
between the behavior of the original Hadamard walk and the unitary noise are clearly shown. 
First, alongside the oscillating wings on the edges, an incoherent component is generated 
around the center [Fig. 12(a)] . Gradually, the central peak increases absorbing weight from 
the two oscillating wings [Fig. 12(b)], which later disappear completely, [Fig. 12(c)] . However, 
the probability distribution does not converge to a Gaussian even at much longer times, since 
at any time step the quantum effects continue to be dominant. The incoherent component 
becomes more dominant as time evolves, due to the accumulation of noise effects. However, 
the coherent features of the quantum walk appear in every time step and the probability 
distribution does not seem to converge to a Gaussian. 

In Fig.Elwe present the probability distribution Pt(n) for the case of strong noise (a = 0.8) 
at t — 1000. At this noise level the quantum effects are completely suppressed by the noise. 
The oscillations are smoothed and the probability distribution converges towards a Gaussian 
distribution with a zero mean, that expands with time. It thus approaches the behavior of 
the classical random walk. 

In Fig. |U we present the time dependence of the standard deviation a a {t) of Ptin) for the 
noisy Hadamard walk with symmetric, unbiased noise. The standard deviation a Q (t), given 
by Eq. (|24j) . is shown for different noise levels a. The top curve is for the noiseless quantum 
walk, the next four curves, from top to bottom, are for a = 0.025, 0.05, 0.1 and 0.2, and 
the last curve at the bottom is for the classical random walk. The initial state used in these 
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simulations is 



\m) = \0)\R)- (26) 

We observe that as the noise level a increases, the standard deviation curves for the noisy 
quantum walk move down from the top curve of the noiseless quantum walk and approach 
the bottom curve describing the classical random walk. 

The standard deviation a Q (t) for the noiseless quantum walk as a function of time t takes 

n 

the form |3|] 

(T (t) = qt. (27) 



From our numerical simulations we obtain that the coefficient q = 0.4505 ± 0.0005, which is 
in perfect „. with the analytical results of Re, 0. The standard deviation for the 
classical random walk is given by 

^classical (*) = ^ ( 28 ) 
As long as the noise level is not too high, the noisy Hadamard walk shows the typical 
behavior of a quantum walk for short times. However, beyond some crossover time the 
motion of the noisy quantum walk acquires diffusive features in the sense that cr a (t) behaves 
like the classical random walk. However, the shape of Pt{n) does not approach a Gaussian 
and maintains a far reaching tail. 

In order to evaluate the crossover time, and examine its dependence on a, we will use 
the short-time and long-time limits of cr a (t). For short times a a (t) increases linearly with 
t, according to Eq. (J27j) . Then, gradually the accumulation of the noise modifies the curve 
until in the long time limit it forms a square-root shape. The long-time tail of u a {t) is well 
fitted by a a (t) = K(a)\/t + C(a), and the functions K(a) and C(a) can be obtained. It 
turns out that for long enough times C(a) can be neglected, giving rise to 

a a {t) = K(a)Vi. (29) 

The function K(a) is then obtained by linear regression of the points [y/i, cr a (t)}, for each 
value of a. The function K(a) for a relevant range of a values is shown in Fig. It 
turns out to be a monotonically decreasing function, that in the limit of small a takes the 
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form K(a) ~ 1/a. In the limit of very large a, K(a) — > 1, and the noisy quantum walk 
coincides with the classical random walk, for which K = 1 (dashed line). The monotonically 
decreasing behavior of K(a) means that as the noise level increases, the broadening of the 
distribution Pt{n) in the long time limit becomes slower. This is due to an earlier suppression 
of the quantum coherence, that broadens the distribution much faster than the noise. 

We will now consider the crossover time T2(a), which is based on the behavior of a a (t). It 
is defined as the time at which the long time asymptote of the noisy Hadamard walk [given 
by Eq. (|29|)] intersects with the short time asymptote, namely the ideal noiseless linear curve 
[Eq. flZZjl ]- This crossover time is given by 

K(a)y/T 2 (a) = qT 2 {a), (30) 

and therefore 

r 2 ( a )=[^l 2 . (3i) 

The crossover time T^cn), obtained from Eq. (j31)l . is shown in Fig. El as a function of a, on 
a log-log scale. It is found that 



T 2 = c 2 a~\ (32) 

where c 2 = 0.62 ± 0.01 and 77 = 2.05 ± 0.08. The crossover time can be interpreted as the 
time required for the noise to scramble the phases and thus eliminate the structure of the 
constructive and destructive interference that characterize the Hadamard walk. This time 
can be estimated by 

av^~l (33) 

namely, 

T 2 ~ i (34) 
cr 

in agreement with the simulation results. Unlike the classical random walk that produces a 
symmetric probability distribution for a particle that starts from the origin, the Hadamard 
walk exhibits an inherent asymmetry. For example, the initial state of Eq. (}2T)|) produces a 
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distribution Pt{n) with much more weight on the right vs. the left hand side. Moreover, the 
average position of the walker moves to the right at a constant speed, namely 



n(i) = vt. (35) 

where v = 0.293 ± 0.003, in perfect agreement with the analytical results of Ref. In 
Fig. [7|we show the evolution in time of the average position (n(t)) a [given by Eq. (|22p] for 
the noisy Hadamard walk with different values of the noise a. The dashed linear line at the 
top is the result for the noiseless Hadamard walk, given by Eq. ([35)1 . Below it, from top to 
bottom are the simulation results for a = 0.025,0.03,0.04,0.07 and 0.1. For short times, 
the average (n(t)) a of the noisy Hadamard walk tends to follow the straight line of Eq. (J35j) . 
up to some crossover time, T\(a), at which it saturates and approaches a constant value: 

(n(t)) a -> n a . (36) 

The asymptotic value, n a , decreases as the noise level a is increased. This is due to the fact 
that at higher noise levels, the buildup of the asymmetric pattern, which is a quantum effect, 
is suppressed more quickly. We define the crossover time Ti(a) as the time t at which the line 
of Eq. (|Hnj), describing the noiseless Hadamard walk, intersects the asymptotic horizontal 
line {n(t)) a = n a , namely vTi = n a . Therefore 

Tx = ^. (37) 

v 

In Fig. Elwe present the crossover time Ti(a) vs. a on a log-log scale. It is well fitted by a 
power-law function of the form 



T-i 



C\ a 



(38) 



where C\ = 0.20 ±0.01 and p = 2.06 ±0.08, which is consistent with the expected asymptotic 
result of p = 2, based on an argument similar to Eq. ()33|) . We thus find that both T\(a) 
and T2(a) exhibit similar dependence on the noise level a, up to a constant factor. 

We will now discuss the connection between our results for the effect of unitary noise on 



the quantum walk and earlier results on decoherence effects 
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23, 



24. 



25, 



Decoherence 



can be described by an additional quantum operation that applies on the chirality qubit at 
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each step of the quantum walk. In a particular choice of the noise model, this operation is 
given by the completely positive map that consists of the operators 

A = ^\R)(R\ 

A 1 = y/p\L)(L\ 



A 2 =y/1- pi, (39) 
where I is the identity operator. This set of operators satisfies 

2 

£4A = /- (40) 

n=0 

It can be interpreted as a measurement of the chirality qubit that is performed with prob- 
ability p. Analytical calculations show that in the long time limit, the first moment of the 
decohered quantum walk, starting in the state |0)|i?), saturates and approaches the value: 



(n(t)) p ^n p = i 77r ^. (41) 



p(2 - p) 

For weak decoherence, or p <C 1, n p ~ 1/p, which resembles our result n a ~ 1/a 2 , if p is 
replaced by a 2 . This is a sensible connection since a represents some matrix elements that 
multiply the amplitudes while p is a probability. 

The results of Ref. for the second moment of the spatial distribution of the decohered 
quantum walk can be expressed by 



a p (t)=K(p)Vi, (42) 

where 



Like K(a), this is a monotonically decreasing function, that in the limit of small p takes 
the form K(p) ~ l/v/p, while K(p) — > 1 for p — > 1. Therefore, associating p with a 2 , one 
obtains the same scaling behavior of the noise effects in the two noise models for both the 
first and second moments of the distribution. 
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V. SUMMARY 



We have studied the effect of unitary noise on the discrete one dimensional quantum 
walk. We have shown that when the quantum walk is exposed to unitary noise it exhibits 
a crossover from the quantum behavior for short times to classical-like behavior for long 
times. For times shorter than the crossover time, the standard deviation a a (t) of the spatial 
distribution of the noisy Hadamard walk, scales linearly with t, like the noiseless Hadamard 
walk. Beyond the crossover time, a a (t) scales like y/t, namely it acquires diffusive behavior, 
like the classical random walk. The crossover time was also characterized using the average 
position of the random walker, namely the first moment, (n(t)) a , of the spatial distribution, 
which scales like t for short times and saturates to a constant value for long times. In both 
cases the crossover time was found to scale as a~ 2 where a is the standard deviation of the 
random noise. 
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FIG. 1: The probability distribution Pt(n) of the noiseless Hadamard walk on the one dimensional 
lattice at time t, vs. the spatial coordinate n. This probability distribution is given for different 
times: t = 250(a), t = 1000 (b) and t = 10000 (c). The initial state of the quantum walk is given 
by Eq. <pg|). 

FIG. 2: The probability distribution (Pt(n)) a of the noisy Hadamard walk at low noise level 
a = 0.025, as a function of n for t = 250 (a), t = 1000 (b) and t = 10000 (c). The initial state of 
the walk is given by Eq. (|25|). 

FIG. 3: The probability distribution {Pt{n)) a of the noisy Hadamard walk for high noise level 
a = 0.8, as a function of n, for t = 1000. The initial state of the walk is given by Eq. (|25l) . 

FIG. 4: The standard deviation a a (t) for the noisy Hadamard walk as a function of the time t, 
for several noise levels a. The top curve (dashed line) is for the noiseless Hadamard walk. Below 
it (from top to bottom) are the results for the noisy Hadamard walk with a = 0.025,0.05,0.1 and 
0.2, respectively (solid lines). The bottom curve (dashed line) is for the classical random walk. 
Each data point was averaged over 200 runs for weak noise (a < 0.07) and over 4000 runs in case 
of strong noise (a > 0.07). 

FIG. 5: The coefficient K(a) as a function of the noise level a. In the limit of small a it takes the 
form K(a) ~ 1/a. In the limit of very large a, K(a) — > 1, namely, it coincides with the result of 
the classical random walk. 

FIG. 6: The crossover time T2 as a function of the noise level a, on a log-log scale. It follows a 
power-law behavior T2 ~ a _r? where 77 = 2.05 ± 0.08. 



FIG. 7: The average position {n(t)) a of noisy Hadamard walk as a function of the time t. The 
top curve is for the noiseless Hadamard walk and below it, from top to bottom are the curves for 
a = 0.025,0.03,0.04,0.07 and 0.1. 

FIG. 8: The crossover time T\ as a function of a, on a log-log scale. It follows a power-law of the 
form T\ ~ a~ p where p = 2.06 ± 0.08 
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